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1. INTRODUCTION 

During the last decade, intensive research on multi-robot systems (MRS) has been conducted to 
solve several real-world problems. Usually, MRS are fault-tolerant and can efficiently solve complex 
problems [1-3]. The issue of task allocation is one of the main challenges in MRS. In such a situation, a 
complex problem is usually divided into elementary and straightforward tasks, where each one is then 
assigned to a robot or several robots [4]. Solving a multi-robot task allocation (MRTA) problem is to assign 
tasks to robots while optimizing a given objective function under several constraints [5]. 

Several centralized (i.e., using a central agent) and decentralized [5-7] approaches have been 
proposed to solve MRTA problems in the literature. Centralized systems reduce information duplication, and 
the check of generated solutions’ quality is effortless [8]. However, maintaining connectivity between robots 
and the central agent is difficult, especially in large real-world problems [5]. Also, a load of communications 
on the central agent can sometimes exceed its tolerable threshold. Besides, centralized approaches are not 
robust (i.e., if the central agent fails, the system fails too). Decentralized approaches generally overcome the 
disadvantages of centralized ones, particularly robustness. However, no guarantee [2] is provided on the 
quality of generated solutions. Currently, we mainly distinguish three taxonomies of MRTA problems, i.e., 
each one classifies MRTA instances according to some features and criteria [9-11]. Taxonomies [9, 10] are 
certainly the most adopted. We discuss them in the next paragraphs. 
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The taxonomy proposed in work [9] counts eight classes: i) [single-task robots, single-robot tasks, 
instantaneous assignment]; ii) [single-task robots, single-robot tasks, time-extended assignment]; 
iii) [single-task robots, multi-robot tasks, instantaneous assignment]; iv) [single-task robots, multi-robot 
tasks, time-extended assignment]; v) [multi-task robots, single-robot tasks, instantaneous assignment]; vi) 
[multi-task robots, single-robot tasks, time-extended assignment]; vii) [multi-task robots, multi-robot tasks, 
instantaneous assignment], and viii) [multi-task robots, multi-robot tasks, time-extended assignment]. In the 
rest of the paper, this taxonomy is called “taxonomy A” for simplicity’s sake. 

—  Single-task robots: This means that each robot can perform a single task at the same time. 

— Multi-task robots: This means that some robots can perform multiple tasks simultaneously. 

—  Single-robot tasks: This means that each task requires exactly one robot for its accomplishment. 

—  Multi-robot tasks: This means that some tasks require the cooperation of several robots for their 
accomplishment. 

— Instantaneous assignment: This means that the available information only allows instantaneous allocations 
of tasks to robots, i.e., we do not use time windows, precedence, or synchronization constraints. 

—  Time-extended assignment: This means that more information is available, i.e., we use time windows, 
precedence, or synchronization constraints, which allows time-extended allocations of tasks to robots; 
in other words, the use of planning is unavoidable. 

The taxonomy proposed in work [10] distinguishes 21 classes, which are: no dependencies (1 or 2), 
in-schedule dependencies (2, 5 or 6), cross-schedule dependencies (1, 2, 3, 4, 5, 6, 7 or 8), and complex 
dependencies (1, 2, 3, 4, 5, 6, 7 or 8). Where the numbers between parentheses represent one of the eight 
classes defined in work [9]. In this taxonomy, each robot r is supposed to have a schedule of assigned tasks. 
In the rest of the paper, this taxonomy is called “taxonomy B” for simplicity’s sake. 

— No dependencies: When allocating a task t to a robot r, we consider t and r. 

—  In-schedule dependencies: When allocating a task t to a robot r, we consider r and its schedule. 

—  Cross-schedule dependencies: When allocating a task t to a robot r, we consider r, its schedule, and 
other robots’ schedules. It is worth pointing out that the robots’ schedules are static, i.e., schedules 
cannot be changed during the allocation process. 

— Complex dependencies: When allocating a task t to a robot r, we considered r, its schedule, and other 
robots’ schedules. It is worth pointing out that the schedules are dynamic, i.e., schedules can be 
changed during the allocation process. 

All MRTA problems are firmly NP-hard [11]. The computational time to generate optimal solutions 
for such problems increases exponentially with the problem size (i.e., numbers of robots and tasks). 
Accordingly, analytical methods are strongly discouraged [12], especially when the problem size is large. In 
the following, we present some related work. 

There exist various common points between the following problems: i) the multi-robot task 
allocation, ii) the multiple traveling repairmen [13], and iii) the multiple traveling salesmen [14]. The last 
problems are derived from the well-known traveling salesman problem [15], which has been demonstrated to 
be NP-hard in combinatorial optimization [16]. Therefore, the multi-robot task allocation problem is also NP- 
hard [5]. The goal of the multiple traveling repairman problem is to determine a tour for each repairer while 
optimizing the average waiting time of customers. The aim of the multiple traveling salesmen problem is to 
allow each salesperson to visit all the cities and then return to the first one while optimizing their traveling 
distances. These criteria (i.e., time and distance) are usually the most adopted in MRTA problems. In the last 
decades, several approaches have been proposed in the literature to solve the different MRTA problems. We 
mainly distinguish between centralized and distributed methods [17-20]. 

In centralized approaches, an intermediate robot collects all the required information and allocates 
tasks to robots. The global objective function is optimized based on complete data [21, 22]. The work 
presented in [23] combines a map partitioning method with a Markovian Decision Process to assign a set of 
exploration tasks to a fleet of robots. Centralized approaches have many advantages: the objective function is 
easily optimized, the simplicity of implementation, and the running time is significantly reduced, to mention 
but a few. Nevertheless, they have some drawbacks, e.g., i) the scalability in real-life applications is not 
always reachable: i.e., guarantee that all robots are connected to the central station leads to a heavy traffic 
burden on it [24]; ii) the coverage area is drastically restricted because the allocations are computed on the 
central station: i.e., the computational demand on the central robot is too high; and iii) the centralized 
approaches are not robust: i.e., single failure point. The work described in [25] proposed a methodology to 
solve a scenario of fire extinguishment that combines quantum genetic algorithms and reinforcement learning. 

In distributed approaches, each robot runs the allocation algorithm individually. Usually, distributed 
methods comprise a consensual phase to guarantee robots’ convergence to consistent allocations, despite the 
network topologies [26]. Consensual stages add supplementary computing charges to the allocation process. 
Consequently, the makespan of these approaches is increased (a significant amount of time is taken by the 
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robots to converge to a consistent allocation) [27]. Auction-based approaches [28-31] are considered as 
distributed methods. They are generally employed to solve problems where each robot can execute one task 
at a time, each task is strictly assigned to one robot, and allocations are time-extended (i.e., ST-SRTA). 
Auction-based algorithms are shown to produce sub-optimal solutions [9] efficiently. In auction-based 
approaches, robots offer bids on tasks according to their respective contextual information. The robot with 
the highest bid wins. The auctioneer might be one of the robots or a dedicated central station [32, 33]. In 
many ST-SR-TA task allocation scenarios, tasks, and robots are tightly bound [34]. In single-round 
combinatorial auctions [35], robots consider task bundles when calculating their bids, and the produced 
allocations are near-optimal. However, the computing time of all possible task bundles is exponentially 
increased with the number of tasks. In sequential single-item auctions [36], this shortcoming is overcome. 
Robots place their bids only on unallocated tasks in each round until allocating all tasks. In sequential 
simultaneous auctions [22], robots place their proposals on unallocated tasks in each round. Only one task is 
assigned at a time to each robot. An unmissable shortcoming in auction-based algorithms is the network 
topology, as the robots must place bids on tasks and transmit them to the auctioneer. The approaches 
proposed in [37-40] presented some auction-based strategies to solve MRTA problems. They used several 
metaheuristics, such as genetic algorithms, firefly algorithm, ant colony system, quantum genetic algorithms, 
artificial bee colony optimization. In work [41], the authors examined the MRTA problem using auctions. 
They analyzed the influence of lossy communications on solutions’ quality. In work [42], the authors 
introduced a method to resolve the MRTA problem. They transformed the situation into a transportation 
problem and proposed a collision-free task allocation scheme. 

In work [43], the authors proposed an algorithm called the “consensus-based bundle algorithm” 
(CBBA). They combined the consensus-based and the auction-based approaches to take their advantages. The 
consensus phase is used to agree on the winning bid values. The algorithm proposed in work [34] is derived from 
CBBA. It is shown to have exceptional performance compared to the CBBA. Other algorithms [44-48] extended 
from CBBA have also been shown to produce better solutions when solving ST-SR-TA problems with 
temporal constraints. This family of algorithms is proven to be very robust, despite the network topology [5]. 
However, they cannot handle the situations where schedules are dynamic. Besides, they produce sub-optimal 
allocations due to greedy-based strategies in the task inclusion phase [5]. The work in [49] uses CBBA, and 
allocations are computed using ant colony optimization. The work described in [50] uses the CBBA, and 
assignments are computed combining ant colony optimization and bat algorithm. 

Many algorithms are available in the literature to solve MRTA problems [51, 52]. The authors of the 
work [53] solved the multi-objective coalition formation problem using the Pareto archived evolution 
strategy algorithm. The proposed method is centralized and does not consider dynamic scenarios. The 
approach described in [18] extended the work in [34] to handle online rescheduling. The algorithm is not 
efficient if much information arrives frequently. The work in [5] presented a solution to the task allocation 
problems in a distributed MRS, under critical time constraints. The authors of the work [54] proposed a multi- 
objective optimization of multi-robot task allocation with precedence constraints. The work in [55] discussed a 
new coordination control for multiple non-holonomic robots for target tracking with limited communications. 

We give a mathematical formulation of the MRTA problem using set theory concepts. It allows a 
formal comprehension of the problem and gives the required tools to develop efficient solutions. At the end 
of the paper, we provide a generic framework for solving any MRTA problem. The work presented in this 
paper has three significant contributions. First, we propose a formal analysis of MRTA problems based on set 
theory concepts. It is not meant to be final or exhaustive, and indeed it has limitations. However, we believe 
that will help researchers understanding the nature of the problem, its components, and its complexity. 
Second, we formalize the taxonomies proposed in [9, 10], using our formal analysis. It will help recognize 
the type of a given problem from its components’ properties (i.e., robots, tasks, allocations.). Third, we 
propose a generic solving scheme (1.e., an algorithm) of MRTA problems. Each step of this scheme is explained, 
and its time complexity is given. Also, we illustrate the proposed generic solving method in a challenging situation, 
which is the assignment of papers to reviewers within a journal. In conclusion, the work presented in this paper is 
not meant to be better than other literature formulations [9-11]; it just complements them. 

The rest of the paper is organized as follows. Some standard definitions specific to set theory are 
given in section 2. Section 3 shows the proposed formal analysis of task allocation in MRS. Sections 4, and 5 provide 
the formulation of different MRTA classes of taxonomies A [9] and B [10], respectively. Section 6 discusses 
the proposed generic solving scheme of MRTA problems and highlights its time complexity. Section 7 
illustrates the proposed algorithm on a real-world problem, which is the assignment of papers to reviewers 
within a journal. Finally, a conclusion and some future work are provided in section 8. 
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2. BACKGROUND ON SET THEORY 

In this section, we give some standard definitions specific to the set theory. 
Definition 1 A set is a collection of objects (i.e., set elements). We write e € E to denote that e is an element 
of E. For example, the set E = {e,, e3, e3} contains three elements. The symbol |E| denotes the cardinality of 
E. The symbol Ø represents the empty set [56] 
Definition 2 A subset is a set whose elements are all members of another set. The term F © E is read: F is a 
subset of E. The term F C E is read: F is a proper subset of E [56]. 
Definition 3 Let E be a set. The power set of E, denoted as (E), is the set that contains all possible subsets 
of E. The cardinality of (E) is 2'E! [56]. The power set of E is given by (1). 


PE) = {X|X S E} (1) 


Definition 4 A collection of subsets F of a set E is called a family of subsets of E. For example, let E = 

{€1,€2,€3,€4,es} be a set, an instance of a family of subsets of E might be given by F = 

ffei, ez, e3}, {€4, es}, {e1,€2,€4}} It is worth pointing out that any family of subsets of E is a subset of the 

power set (E) if it has no repeated members [56]. The cardinality of a family of subsets F is equal to the 

number of subsets it contains. The terms “family of subsets F ” and “family F ” will be used interchangeably. 

Example 1 Let E = {e}, ez, e3} be a set. 

— The power set of E is: (E) = {ø, {e1} {22} {23} fer ez}, fer, es}, {e2, es}, {er e2, e3}}- 

— Collections of subsets F} = ffe, ez}, {€2, ez3}}, F, = {fei}, {e} and F}; = ffe, ez, es}} are three families 
of E. The cardinality of F, is |F2| = 2. 

Definition 5 The Cartesian product E; X ...X En of sets E1, ..., En is defined as the set of all ordered tuples 

(e1, =, €n), with ey E Ey,...,€, E En as follows [56]. 


E; X X En = {(ey, en) le, E Ey, €n E En} (2) 


Definition 6 Let E}, ..., En be n sets. An n-ary relation on these sets is a subset of E, X ... X E,. The sets 
E,,...,E, are called domains of the relation, and n is its degree [56]. 

Definition 7 A k-permutation is an ordered arrangement of k subsets from a family of n subsets (k < n). Ifa 
given family of subsets F has n subsets, then the number of k-permutations generated from F, denoted as 


P(n,k), is given by (3) [56]. 


P(n,k) = a (3) 


Example 2 Let F = {fe}, {ez}, {e3}, {e1, en}} be a family. It is evident that F has 24 different 4-permutations, 
e.g., ffei}, {ez}, {e1, eo}, {e3} is one of them. 

Definition 8 Let E' be a subset of a set E and F' be a subset of a set F. According to work discussed in [10], 
we can define a utility measure, denoted as U (E', F"), for subsets E' and F', (4). 


V(E',F') , if elements of E'and F'are related 
00 , otherwise 


UCE", F’) = { (4) 


The term V(E',F') is any function that returns a numerical value according to relationships between 
elements of subsets E' and F'. For example, if elements of E' and F' are points in the space, V(E', F") may 
return the average of distances between points of E' and points of F'. In other words, the function V (E', F") is 
given (5). 


E F 
LEI eles ll 
|E"|.[Fr| 


V(E',F') = (5) 


3. FORMAL ANALYSIS OF MRTA PROBLEMS 

In this section, we give the proposed formal analysis of multi-robot task allocation problems. 
Definition 9 A robot is an autonomous (logical or physical) entity that perceives its environment, makes 
decisions, and performs actions [11]. The set A = {Q4, ..., An} designates a group of n robots. 
Definition 10 A task is an elementary and independent action accomplished by one or several robots [11]. 
The set T = {ty, ..., tm} designates a collection of m tasks. 
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Definition 11 A resource refers to any tool (logical or physical) that a task requires for its accomplishment. 
For example, suppose we want to measure the ambient temperature of someplace (i.e., task). In that case, we 
need a thermometer (i.e., resource). In MRTA scenarios, robots usually offer resources and tasks use them. 
The set Q = {@y, ..., Wg} designates a collection of k resources. 

Definition 12 Resources offered by robots is a ternary relation over sets A, Q, and B = {0,1}. It is a set of 
ordered triples (a,w,b) consisting of elements: a€ A, w EQ, and b EB. It encodes the relation’s 
information: robot a has b instances of resource w. This relation is denoted as O (O & A x Q x B). 
Definition 13 Resources required by tasks is a ternary relation over sets T, Q, and B = {0,1}. It is a set of 
ordered triples (a,w,b) consisting of elements: tET, w EQ, and b EB. It encodes the relation’s 
information: task t requires b instances of resource w for its accomplishment. This relation is denoted as R 
(RETXO2XB). 

Definition 14 Valid assignments over robots and tasks is a quinary relation over sets A, B = {0,1}, 0, T, and 
B = {0,1}. It is a set of ordered quintuples (a, b, w, t, b') consisting of elements: a EA, b €B,w EQ, t ET, 
and b' E€ B. It encodes the relation: robot a has b instances of resource w and can be assigned to task t, 
which needs b' instances of resource w. This relation is denoted as Y (VG AXBxQxT Xx B). 

This relation is given by (6), where Q;(a, b, w,t, b') are conditions that determine allowed quintuples. For 
example, suppose the distance separating a robot and a task is larger than a given threshold. In that case, 
this robot cannot be assigned to this task even if it has the required resources. In other words, these 
conditions translate the constraints of the considered problem. 


V = {(a,b,, t, b')|(a, w, b) E O, (t, w,b') E R, (b,b’) = 
(1,1), Q, (a,b, w, t, b'), Q2(a, b, w, t, b’), Q3(a, b, œ, t, b’), ...} (6) 


NB: The set B = {0,1} can be replaced by N = {0,1,2,3, ... }. However, for simplicity, we use the first set. 
Definition 15 A coalition is a temporary group of robots working together on some task. Generally, when the 
task is performed, the coalition is dissolved [11]. If we consider a set E, then one can enumerate all possible 
coalitions using Definition 3. In this context, we should exclude the empty set ® since it is not a coalition. 
Thus, for a given set E, we have 2'*! — 1 possible coalitions. 

Definition 16 Let E be a set. An overlapping family F of E is a family whose subsets are not duplicated and 
not necessarily pairwise disjoint. An overlapping family F should simultaneously verify the following 
conditions: 


a Ø ¢F. 
b. aF, F, € F|(F, # Fy) > (F, N F, + Ø). 
c. UË F, =E. 


Property 1 signifies that the empty subset is not considered (i.e., the empty coalition is not considered). 
Property 2 implies that a given element can be a member of several subsets simultaneously (i.e., a given 
robot can belong to several coalitions at the same time). Property 3 denotes that all members of set E should 
appear in subsets of family F (i.e., all robots are allocated). The symbol O(E) designates the collection of all 
possible overlapping families of subsets of a set E, the symbol ©; denotes the i" family in the collection @(E), 
and the symbol ©;; expresses the j subset in the family ©;. The number of overlapping families of a set E is 
given by (7) [57]. 


|El! 2lEl-k 


IO(E)| = han Go (|E|—k)!k! 





(7) 


Definition 17 Let E be a set. A disjoint family F of E is a family whose subsets are not duplicated and 
mandatory pairwise disjoint. A disjoint family F should simultaneously verify the following conditions: 

a @O€ F. 

b. VFF, E F|(F + F) > (Fy NF, = Ø). 

c UEF; =E. 

Property 1 signifies that the empty subset is not considered (e.g., the empty coalition is not considered). 
Property 2 implies that a given element cannot be a member of several subsets simultaneously (e.g., any task 
is allocated one time). Property 3 denotes that all members of set E should appear in subsets of family F 
(e.g., all tasks are allocated). The symbol A(E) designates the collection of all possible disjoint families of 
subsets of a set E, the symbol A; denotes the i family in the collection ACE), and the symbol Aij expresses 
the j" subset in the family A;. The number of disjoint families of a set E is given by (8) [58]. 


k! 
(k-i)!i! 





IME) = Dito (FEO gg D) i 
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Remark 1 Jf for each family A,;€ ACE) we consider its \A;l-permutations, we can define an augmented 
collection A' (E), where each family A;€ A(E) is replaced by its \A,\-permutations. The symbol A;' will 
denote the i" family in the collection A'(E) and the symbol Ajj’ will express the ji" subset in the family Aj‘. 
The number of augmented disjoint families of a set E is given by (9). 





|A’(E)| = Deco (G OD a 0) r) (9) 


Example 3 Let E = {e}, e2} be a set. 
— The collection of all possible overlapping families is: 


OCE) = f{fe,}, {e2}, {Ler e2}, {fer}, fer ead} {feck fer, ex}}, {fer}, fer} fer, ex} }}. 
— The collection of all possible disjoint families is: 

ACE) = {{fe1}, fez}, {fer e}}}. 
— The collection of all possible augmented disjoint families is: 


A'(E) = {{fe.}, fe2}}, {fea}, {e13}, {fer e23}}. 


Definition 18 Let A = {a,,...,a,} be a set of robots and T = {t,,...,tm} be a set of tasks. Let O(A) = 
{0,, On} be the collection of overlapping families of A and A'(T) = {A}, mo Ag} be the augmented 
collection of disjoint families of T. Terms O; and A,’ are, respectively, called family of robots’ coalitions and 
family of task’ coalitions. 

Partial allocations are a binary relation over collections @(A) and A'(T). It is a set of ordered pairs 
(;, A;’) consisting of elements: ©; E€ O(A) and A,’ € A'(T). It is a subset of the Cartesian product O(A) x 
A'(T). This relation is given by (10) and denoted as T. 


T = {(0;,A,')|@; € O(A), Aj’ € A’ (T), 10; = |A;"|, Q1 (03), «., Q1'(Aj’), s Qa (O Aj"), (10) 


Each pair (0;,A,') ET is read as follows: each robots’ coalition Oig E ©; is assigned to tasks’ coalition 
Ajk € Aj’. For example, let 0; = {{a,, a3}, {a a2} and Aj= fit, t2}, {t3}}. We read: robots’ coalition 
{a,, a3} is assigned to tasks’ coalition {t,,t}, and robots’ coalition {a}, Az} is assigned to tasks’ coalition 
{tz}. Q (00), =, Q,'(A;’), en Q,''(0;,4,'), ... are conditions that determine allowed partial allocations. 
Definition 19 A global allocation T' (or merely an allocation) is a sub-relation of T whose pairs (®;,A;') 
should simultaneously verify the following conditions: 


UM Ulo, =A. 
b O UR Ay’ ST, 

AI aR Ay’ =Ø. 
d. V(@;,A;') E T'YO;j E OVALE Aj|{wlw E (O; x Bx Ax Aj; x B)} N {wlw E (Aj; x 2x B)} = 

{wlw E (Aj; x Ox B)} 
Relation (Ojj xXBXxXAX Aij x B) E V represents valid assignments over robots in @;; and tasks in Aij- 
Relation (Ai; xQ x B) E R represents the required resources by tasks in A;j. Condition 1 denotes that all 
robots are solicited to be assigned to tasks. Condition 2 means that all tasks should be assigned. Condition 3 
signifies that a given task is exactly assigned to one robots’ coalition. Condition 4 states that each coalition 
of robots should offer the required resources to the corresponding coalition of tasks. The optimal global 
allocation is the one that optimizes the utility measure given by (11). 

U= ll u(y, Ay’ (11) 
Definition 20 Let (F,, F2) ET be a partial allocation. Let F}; and F3; be the i robots’ and tasks’ coalitions 
of families F, and F,, respectively. Potential assignments are a sub-relation of the relation V given by 
Definition 14 (valid assignments) whose quintuples modeled by (12). This relation is denoted as P. 


P = {(a,b,a,t, b’)|(a,b,,t, b') E V,a E Fii t E Fay} (12) 
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Robots in F,; ar related to tasks in Fz; (4), if and only if the relation P is not empty. 

Example 4 In this example, we consider A = {a,,a2}, T = {t1, t2, t3}, and Q = {w,, wz}. Relations O, R, 
and Y are given as follow: 

- O= {(4,, w, 1), (a, %2,0), (a2, w, 1), (a2, w2, 1)}. 

a e az {t Wi, 1), (ti, w2, 0), (tz W1, 1), (tz, W2, 0), (tz, Wi, 0), (ts, W2, 1)}. 

= = = {(a, 1, W1, ti, 1), (a1, 1, W1, t2, 1), (az, 1, W1, ti, 1), (az, 1, W1, tz, 1), (az, 1, W2, tz, 1)}. 

Let = {fa az}, {a}} be a family of robots’ coalitions; and F, = {ft}, {tz t3} be a family of tasks’ 
coalitions. If we consider F} = {az} and Fy, = {tz, tz}, then P = {(a,, 1, w1, tz, 1), (az, 1, w2, tz, 1)}. 
Definition 21 Let T’ = {(E,, F4), ..., (Ep, Fp)} be a global allocation. A family of potential assignments is the 
set [P] = {P,, Ley Pp} of potential assignments given by Definition 20 (each term P; € [P] is computed as in 
Example 4). 

The concept of a family of potential assignments allows us to build task schedules of robots. The use of this 
concept will be detailed in Section 5. 


4. FORMULATION OF TAXONOMY A 

This section gives the formulation of taxonomy A [9], using our proposed formulation. 
Definition 22 Let © be the relation of resources offered by robots given in Definition 12. All robots of a 
given problem are said to be Single-Task, if the following condition is satisfied. 


Va € AJ! w E Al(a,w, 1) € O (13) 


Condition 13 signifies that each robot a € A offers one and only one resource w E Q for accomplishing tasks. 
Definition 23 Let © be the relation of resources offered by robots given in Definition 12. Some robots of a 
given problem are said to be Multi-Task, if the following condition is satisfied. 


Va € Ady, ..., 0p E A\(a,@1,1) EOA...A (a, @py,1) € O (14) 


Condition 14 signifies that some robots a E A offer several resources w E Q for accomplishing tasks. 
Definition 24 Let R be the relation of resources required by tasks given in Definition 13. All tasks of a given 
problem are said to be Single-Robot, if the following condition is satisfied. 


Vt ETAlw E A|(t,w,1) E R (15) 


Condition 15 signifies that each task t ET requires one and only one resource w EQ for its 
accomplishment. 

Definition 25 Let R be the relation of resources required by tasks given in Definition 13. Some tasks of a 
given problem are said to be Multi-Robot, if the following condition is satisfied. 


Vt E TAQ, ...,@p E Al(t,@1,1) E RA A (t, @p,1) E R (16) 


Condition 16 means that some tasks t E T require several resources w E Q. for their accomplishment. 

Definition 26 Let I’ = {(F,, F,’), sie (F, LE n} be an allocation. 
If the execution of tasks in rv are associated with any temporal model (e.g., time intervals [59], simple 
temporal networks [60]). In this case, we say that the allocation of robots to tasks is performed in a 
time-extended manner. 

— Conversely, we say that the allocation of robots to tasks is performed in an instantaneous manner (i.e., 
no temporal model is provided, and tasks are allocated as they come).’ 


5. FORMULATION OF TAXONOMY B 

In this section, we formulate the taxonomy described in [10], using the proposed analysis. It is worth 
pointing out that each robot should have a schedule of tasks, which is a critical issue in this context. The 
following example shows how robots’ schedules are built. 
Example 5 In this example, we consider A = {a,, a2, a3}, T = {t,, tz, t3, t4}, and Q = {w1, wz}. Relations O, 
R, V, and I" are given as follow: 


ia O= {(a1, W1, 1), (a1, W2, 1), (az, W1, 1), (az w2, 1), (az, W1, 1), (az, W2, 1)}. 
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ki = be w1, 1), (ty, w2, 1), (tz, w1, 1), (to, W2, ma 


(ts, W1, 1), (tz, w2, 1), (ty, W1, 1), (ta, W2, 1) 
(az, 1, W1, ty, 1), (ay, 1, W1, to, 1), (a; 1, W1, tz, 1), (a1, 1, W1, ty, 1), 


(a1, 1, W2, ti 1), (a1, 1, w2, tz, 1), (az, 1, W2, tz, 1), (a1, 1, W2, t4, 1), 

(az, 1, w1, ti 1), (az, 1, w1, to, 1), (az, 1, w1, tz, 1), (az, 1, w1, ty, 1), 

(az, 1, Wo, ty, 1), (az, 1, %2, t2, 1), (az, 1, Wo, tz, 1), (az, 1, Wo, ty, 1), [7 

(az, 1, %1, t,, 1), (a3, 1, 04, tz, 1), (az, 1, w, tz, 1), (az, 1, %1, ty, 1), 

(az, 1, W2, ty, 1), (az, 1, We, tz, 1), (az, 1, We, tz, 1), (az, 1, We, t4, 1) 

- f= (at, p (lta, Az}, {a1, Q2, az}}, {{t2}, {ts, t4}})}. 

T” is a global allocation, which means that: robots’ coalition {a,} is assigned to tasks’ coalition {t,}, robots’ 
coalition {a}, Az} is assigned to tasks’ coalition {t3}, and robots’ coalition {a1, az, a3} is assigned to tasks’ 
coalition {t3, t4}. The family of potential assignments corresponding to T" is [P] = {P4, P2, P3}, where: 

— Potential assignments P, over subsets {a,} and {t,} are: 


- V= 


P, = {(ay, 1, w1 t1, 1), (ai, 1, W2, ty, 1)} 


Robots’ coalition {a,} is assigned to tasks’ coalition {t,}. We have only one partial schedule: i.e., robot a, 
offers both resources w and w; to task ty. 
— Potential assignments P, over subsets {a,, a} and {t3} are: 


P2 = {(a, 1, W1, tz, 1), (a1, 1, W2, tz, 1), (az, 1, W1, tz, 1), (az 1, W2, t2, 1)} 


Robots’ coalition {a}, az} is assigned to tasks’ coalition {t3}. We have four partial schedules: i.e., robot a, 
offers both resources œw, and w, to task tz; robot az offers both resources w and Wy to task tz; robot a, 
offers resource 0), to task tz and robot a, offers resource œw, to task tz; or robot a, offers resource œw, to 
task tz and robot az offers resource w; to task tz. 

— Potential assignments P, over subsets {a1, Az, a3} and {t3,t,} are: 


(a, 1, 4, t3, 1), (ai, 1, w72, t3, 1), (ay, 1, 4, ty, 1), (ay, 1; w72, ty, 1), 
P, = 4 (a2, 1, Wy, t3, 1), (az, 1, W2, t3, 1), (az, 1, %1, ty, 1), (az, 1, We, t4, 1), 
(a3, 1, w1, t3, 1), (a3, 1, w2, tz, 1), (a3, 1, w1, ty, 1), (a3, 1, w32, t4, 1) 


Robots’ coalition {a,, a2, a3} is assigned to tasks’ coalition {t3, t4}. We have 3?%? = 81 partial schedules 


(we only cite one instance): e.g., robot a, offers resource œ to task t}, robot a, offers resource œ, to task 
t3, robot az offers resource œ to task t4, and robot a3 offers resource w, to task t4. 
Generally, if we have three sets A, T and Q. of robots, tasks, and resources, respectively. In this case, we 
should consider \|A|""'*!°D different partial schedules of size \T| x |Q. In our example, we should consider 
1 x 4x 81 global schedules, and the schedule of each robot will contain its assigned tasks: e.g., if we take 
the partial schedules: i) robot a, offers both resources œw and œ; to task t,; ii) robot aq offers resource 04 
to task t, and robot a, offers resource œ, to task tz; and iii) robot a, offers resource œw; to task t3, robot a, 
offers resource w, to task t3, robot a, offers resource œ; to task t4 and robot a; offers resource œ, to task 
t4. In this case, i) the schedule of robot a, contains tasks {(t,, w1), (ty, w2), (tz, w1), (tz, w1), (t3, W2)}; ii) 
the schedule of robot a, contains tasks {(tz, w2), (t4, @1)}; and iii) the schedule of robot a3 contains one 
task {(t4, @2)}. It is worth pointing out that the order of execution of tasks of each robot depends on how 
tasks are allocated to robots, i.e., instantaneous or time-extended assignments. 
Definition 27 Let T' = {(E1, Fa), SAG (Ey, Fy)} be an allocation and [P] = {P,, Jk Pp} be its corresponding 
family of potential assignments. We consider robots a, ...,an. We showed, via Example 5, how robots’ 
schedules are built. Let Sq,(t) be the schedule of robot a; at time t. We define the features of taxonomy B 
[10] as follows. 
— No dependencies: 

As shown in (17) shows that a given robot’s utility measure for a specific task depends on them. 


Vt; E Sa; (0): U({a;}, {t;}) = fla tj) (17) 


—  In-schedule dependencies: 
As shown in (18) shows that a given robot’s utility measure for a specific task depends on its schedule. 


Vt; € Sq,(0): U({a;}, {t;}) = f (ai, Sa, (0)) (18) 


— Cross-schedule dependencies: 
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As shown in (19) shows that a given robot’s utility measure for a specific task depends on all robots’ 
schedules, including itself. However, the schedules are static. 


Vt; € Sq,(0): U({a;}, {t;}) = f (a;,52,(0), TO) (19) 


— Complex dependencies: 
As shown in (20) shows that a given robot’s utility measure for a specific task depends on all robots’ 
schedules, including itself. However, the schedules are dynamic. 


Vt; € Sa; (0): U({ai},{t)}) = f (ai Sa, (T), anes San (x)) (20) 


6. SOLVING SCHEME OF MRTA PROBLEMS 
This section overviews how the proposed formulation is used in a generic solving scheme and 
explains its analytical time complexity. Algorithm 1 highlights the main steps in solving a MRTA problem. 

— Line 1 computes the relation V, a subset of the Cartesian product of relations © and R. The relation V 
contains all valid assignments between tasks and robots. The time complexity of this step is 
O(max{|O|, |RI}*). 

— Line 2 computes the collection O(A) of overlapping families of subsets of the set A; i.e., a given robot 
may belong to several coalitions, since it can perform several tasks (simultaneously or sequentially). The 
time complexity of this step is O(I@(A)| x Ziel), 

— Line 3 computes the collection A(T) of disjoint families of subsets of the set T; i.e., a given task must 
belong to only one coalition, since it must be done once. The time complexity of this step is also 
O(l@(T)| x 2le@l), 

— Line 4 computes the collection A’(T) of all possible permutations of families in A(T). The time 


: ; : Äi 
complexity of this step is O ( max (lA; | n). 


— Line 5 computes the relation I’, a subset of the Cartesian product @(A) x A'(T). This relation contains all 
partial allocations. The time complexity of this step is O(max{|@(A)|, |A’(T)]|}*). 
— Line 6 computes all global allocations. The time complexity of this step is O(ler )| x 2! I, 


— Finally, each global allocation is evaluated using the considered utility measure subject to some 
constraints. The optimal assignment of tasks in T to robots in A is returned. 


Agorithm 1. Generic solving of a MRTA problem 


Input: Gn}, i.e. set of robots, given by Definition 9 
Input: tm}, i.e. set of tasks, given by Definition 10 
Input: wr}, ie. set of resources, given by Definition 11 
Input: O, i.e. resources offered by robots, given by Definition 12 
Input: R, i.e. resources required by tasks, given by Definition 13 
Input: U, ie. utility measure, given by Definition 8 

Input: {<),...,€,}, i.e. constraints to be satisfied 


Output: Optimal assignment of tasks in T to robots in A 


1 Compute V, i.e. valid assignments over robots and tasks, given by Equation 6; 
2 Compute O(A), i.e. the collection of overlapping families of subsets of A, given by Definition 16; 
3 Compute A(T), i.e. the collection of disjoint families of subsets of T, given by Definition 17; 
4 Compute A’(T), i.e. the collection of augmented disjoint families of subsets of T, given by Remark 1; 
5 Compute I, i.e. partial allocations, given by Equation 10; 
6 Compute P(T), i.e. the collection of subsets of the relation T; 
7 foreach I’ € P(T) do 
8 if I’ is a global allocation, given by Definition 19 then 
Compute [P], given by Definition 21; 


Evaluate I’ using U, given by Equation 11; 
Save the optimal assignment; 
end 





It is worth pointing out that MRTA problems are strongly NP-hard, i.e., shown through complexity 
analysis. Therefore, exact methods are not highly recommended to handle large-scale MRTA problems. 
Thus, the use of metaheuristics-based techniques is recommended [37, 49]. 
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7. ILLUSTRATIVE EXAMPLE 

We give an illustrative example to show how Algorithm 1 is employed to solve a given MRTA 
problem. We choose the problem of assigning papers to reviewers within a journal. We suppose that the 
Bulletin of Electrical Engineering and Informatics has received 20 papers to consider for publication on a 
given day. The editor-in-chief (EIC) has ten potential reviewers. After checking all the papers, the EIC 
concluded that they are articulated around five topics. Hence, he should assign each paper to some reviewers 
(e.g., two reviewers) while optimizing a given criterion under some constraints. This situation is indeed a 
MRTA problem, where: (a) reviewers are the robots (Definition 9), i.e., |A| = 10; (b) papers are the tasks 
(Definition 10), i.e., |T] = 20; and (c) topics are the resources (Definition 11), i.e., |Q] = 5. Topics for which 
reviewers are experts are summarized in Table | (Definition 12), i.e., relation ©. Topics covered by papers 
are summarized in Table 2 (Definition 13), i.e., relation R. 

Line 1 of Algorithm 1: We compute V, i.e., valid assignments over robots and tasks, using (6) and 
Tables 1 and 2 (Definition 14). We suppose that no constraint is considered: i.e., {é,, peer) is empty (all 
valid assignments are accepted). It is worth pointing out that in the worst case, we should have |A| x |T| x 
|Q] valid assignments: i.e., |V| = |A| x |T| x |Q]. For our example, |V| = 10 x 20 x 5 = 1000, which is a 
large number to consider. 

Line 2, 3, and 4 of Algorithm 1: We compute @(A) (Definition 16), A(T) (Definition 17), and A’(T) 
(Remark 1), respectively. Cardinalities of collections @(A), A(T), and A’(T) are given as follows. There is a 
tip to reduce the size of collections @(A) and A’(T). For instance, if the journal policy states that each paper 
is assigned exactly to 2 reviewers; therefore, the EIC might consider: (a) families ©; E O(A) for which 
|©; | = 2; and (b) families A,’ € A’(T) for which |A;;'| = 1. 


i 10! 10-k 
|OCA)| = 52k CD" Gop 


IAI = D220 (EEOC Di  (k - 0°) 


(k-i)!i! 


a'r) = AE MCU gage D70) r) 


(k-i)!i! 


Line 5 and 6 of Algorithm 1: We compute I = @(A) x A’(T), partial allocations, given by (10) and 
(I), the powerset of relation T. Line 7 to 13 of Algorithm 1: For each subset in so(1), we check if it is a 
global allocation (Definition 19). If this is the case, we evaluate the global allocations and return the optimal 
one. For the current example, line 9 is ignored since assignments are instantaneous. For example, if we have: 


Aae a7}, {t1}), (fao, a7}, {to}), Cas a10} {ts}), (fag a1} {tah), Caz, a10} {ts}. 
r= ({ag, a3}, {te}), (fag, as} {t7}, Qaz ao), {ts}, (az, a10) {to}), (ar, a4} {t10}, 
({a4, ao}, {t11}, as, ao}, {t12}, az ach {tish), (fae, as} {t14 (Aac a2} {tis D, 
({ag, a7}, {t16), (fas, aro}, {t17), (lao, aro}, {tigh), Caz, ae}, {t19}), Claro, ao} {t20}) 


We evaluate I’ using the objective function given by (21). We sum portions of topics offered by two 
reviewers compared to topics needed by a paper. The optimal global assignment is the one that maximizes 
this sum. Of course, the EIC might use other criteria to evaluate a global allocation: e.g., impose a maximum 
occurrence of a given reviewer. 





g n yl lait [ee (21) 


aa |{wlwe(ajjxoxB)}] 


Table 1. Topics of reviewers (relation O) 








ay Q2 Q3 Qy Qs Qe Q7 Qg (Qo Qio 
w1 1 1 0 1 0 1 1 0 1 1 
wz 1 0 0 0 1 0 1 1 0 0 
w 0 0 0 1 0 1 1 0 1 1 
wo, 0 0 1 1 1 1 0 0 0 0 
Ws 1 1 0 0 1 1 0 0 0 1 
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Table 2. Topics covered by papers (relation R) 








t t ts t ts te t7 tẹ tọ to tu t2 t3 tọ tis tie tı7 tig tio too 
w1 1 1 0 1 0 0 1 0 1 0 1 1 1 1 0 1 1 0 1 1 
wz 1 1 0 1 1 0 1 0 1 1 1 0 1 1 0 0 0 0 1 0 
w, 0 1 0 0 0 1 0 0 1 1 0 1 1 1 1 0 1 0 1 1 
w 0 1 1 0 0 0 1 1 0 1 1 0 1 1 1 1 0 0 1 0 
ws, 0 1 0 1 0 0 0 0 1 0 0 0 0 0 0 1 1 1 0 0 





8. CONCLUSION AND FUTURE WORK 

This paper proposed a formal analysis of multi-robot task allocation problems using elementary set 
theory concepts. We formulated two well-known taxonomies of MRTA problems (i.e., the taxonomy of 
Gerkey and Matari¢ entitled “A formal analysis and taxonomy of task allocation in multi-robot systems”; and 
the taxonomy of Korsah, Stentz, and Dias entitled “A comprehensive taxonomy for multi-robot task 
allocation”), then we developed them using the proposed analysis. A generic solving scheme of MRTA 
problems is presented, and its time complexity is discussed. Finally, the proposed generic solving method is 
employed to solve a real-world problem: i.e., the assignment of papers to reviewers within a journal. In 
conclusion, the work introduced in this paper is not intended to be better than the previous formulations; it 
merely complements them. As future work, we plan to represent the taxonomy of Ernesto, Marie, Hakim, and 
Maria entitled “A taxonomy for task allocation problems with temporal and ordering constraints” using our 
framework. 
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